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INTRODUCTION 

Background 

In  the  course  of  many  physical  measurements,  a  relatively  small  collection 
of  raw  data  is  used  in  a  series  of  calculations,  transformations,  graphs  and 
data  reduction  routines  until  the  final  result  consisting  of  a  parameterized 
equation,  a  graph  or  a  table  is  produced.  In  the  course  of  the  calculations 
it  is  often  very  desireable  to  use  the  same  operation,  or  series  of 
operations,  on  each  set  of  observations  and  similar  calculations  are  often 
performed  on  unrelated  types  of  data.  A  careful  investigator  often  wants  to 
list  and  check  intermediate  data,  or  perhaps  examine  it  using  some  testing 
algorithm,  perform  some  systematic  corrections,  and  keep  the  data  in  some  safe 
and  convenient  form  until  his  results  have  been  confirmed  by  the  scrutiny  of 
others. 

Before  the  advent  of  the  computer,  many  of  the  calculations  involving 
minor  corrections  or  repetitive  measurements  were  simply  not  done.  A  minimal 
amount  of  data  was  taken  and  processed.  There  were  frequently  questions  as  to 
whether  all  observations  had  been  treated  consistently  during  the  reduction 
process,  and  much  time  and  effort  was  often  spent  in  investigating 
discrepancies  that  resulted  from  calculaJ:ional  error  during  data  reduction. 
Even  with  the  aid  of  the  computer,  calculations  were  often  done  piecemeal  in  a 
series  of  programs,  or  a  special  program  to  perform  a  series  of  calculations 
was  laboriously  built  for  each  individual  experiment. 

The  concept  of  a  database  management  system  has  allowed  the  investigators 
in  this  laWKlfiifory'to  collect  and  process  data  systematically,  and  with 
greater  confidence  in  calculations!  consistency,  since  early  1977.  This 
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database  management  system,  called  PHD  (Program  for  Handling  Data),  has  been 
developed  to  support  the  collection  of  physical  measurement  data.  This  type 
of  data  usually  consists  of  a  small  number  of  observations  (2-100),  each 
observation  consisting  of  a  set  of  measured  variables  and  their  uncertainties 
in  some  specified  set  of  units.  For  convenience,  an  observation  identifier 
may  be  included.  PHD  was  developed  to  support  this  type  of  data  collection. 

PHD  consists  of  a  series  of  programs  written  primarily  in  BASIC  which 
allow  the  manipulation  of  data  mathematically.  In  addition,  it  provides 
visual  aids  such  as  listings  and  graphics  as  well  as  the  facility  to  add, 
delete,  correct,  sort,  merge,  rearrange,  store,  retrieve  and  perform 
multivariate  non-linear  regression  on  collections  of  data  stored  as  files 
under  the  host  computer  operating  system  (RSXll-M  on  the  PDP  11/45  or  11/10, 
MINC  BASIC  on  the  MINC,  and  VAX/VMS  on  the  VAX  11/780).  The  system  is 
designed  to  allow  users  to  add  their  ov/n  manipulation  programs  to  those 
available  in  the  system.  The  following  document  describes  the  capabilities  of 
the  software  available  in  the  system  and  how  to  use  them;  it  also  provides 
sufficient  detail  to  allow  users  to  write  routines  for  additional  processing 
of  data  while  preserving  the  advantages  of  the  facilities  available  to  manage 
data  sets  within  the  system.  While  some  of  the  programs  are  peculiar  to  the 
type  of  hardware  being  used — the  graphics  routines,  for  example — the 
algorithms  used  are  adaptable  to  many  systems  with,  as  the  saying  goes,  only 
minor  changes  in  the  coding. 

PHD  is  a  modular  system;  the  main  program  defines  the  internal  data 
storage  system  and  then  passes  control  to  a  selector  or  menu  program,  MENU. 
MENU  calls  each  module  requested,  including  user  written  programs,  and  each  of 
the  modules  it  calls  normally  exits  back  to  the  MENU.  The  description  of  the 
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data  contained  within  the  database  (called  the  BASE)  is  maintained  in  main 
memory  along  with  some  ancillary  storage  used  in  data  manipulation.  The  data 
in  the  database  are  maintained  in  either  main  memory  or  a  virtual  array  which 
is  stored  primarily  on  a  random  access  mass  storage  device  (e.g.,  a  disk). 
When  necessary  the  MENU  cleans  up  the  disk  by  deleting  the  virtual  array  or 
the  temporary  transfer  files  from  the  mass  storage  device  when  they  are  no 
longer  needed. 

Because  of  the  inability  to  trap  errors  in  some  versions  of  BASIC, 
occurrence  of  errors  often  results  in  what  can  easily  be  a  catastrophic  loss 
of  the  data  currently  being  used.  The  system  has  been  designed  to  compensate 
for  this  by  exposing  only  the  data  in  the  current  WORKSET  to  loss  from  an 
untrappable  error,  and  generally  recovery  of  data  after  such  an  error  is  a 
simple  process. 

Acknowledgements 

A  number  of  people  have  contributed  to  the  conception  and  testing  of  this 
program  and  its  various  modules.  I  would  like  to  express  my  thanks  to  all. 
Deserving  special  thanks  for  their  involvement  in  this  process  are  Dr  Lowell 
A.  King  and  Dr  John  S.  Wilkes.  Lt  Robert  E.  Wheelock  concieved  the  algorithm 
and  coded  the  initial  version  of  the  plotting  module  in  1978.  The  bulk  of 
this  module  remains  intact  in  the  current  version.  1  also  would  like  to  thank 
Lt  Col  Chester  J.  Dymek  for  very  helpful  discussions  and  suggestions  on  the 
docustentat  ion. 
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Concept 

Each  set  of  data  (DATASET)  in  the  PHD  system  may  be  thought  of  as  a 
two-dimensional  array  of  data  with  each  observation  occupying  a  row. 

Therefore  the  number  of  rows  in  the  array  is  equal  to  the  number  of 
observations.  Associated  with  each  observation  (or  with  each  row  in  the 
array)  is  an  observation  identifier,  composed  of  two  parts,  a  plot  symbol 
number  (PSN)  and  an  alphanumeric  identifier  (ID).  The  columns  of  the  array 
correspond  to  an  ordered  set  of  variables  measured  for  each  observation. 

Each  column  has  associated  with  it  a  NAME  for  the  variable,  the  UNIT  in 
which  the  data  is  measured,  and  a  value  related  to  the  uncertainty  with  which 
this  variable  is  measured,  the  METRIC.  In  addition  each  DATASET  has  a 
one-line  COMMENT  stored  with  it.  The  first  part  of  this  COMMENT  is  generated 
by  the  system  and  consists  of  the  date  and  time  at  which  the  DATASET  was  SAVEd 
(i.e.,  made  a  permanent  file  on  the  mass  storage  device).  After  it  is  SAVEd, 
a  DATASET  may  not  itself  be  modified.  As  a  safeguard  against  accidental  loss 
of  data,  it  resides  unmodified  on  the  mass  storage  device  until  explicitly 
deleted.  Only  the  current  or  working  dataset  (WORKSET),  which  exists  only 
during  the  active  operation  of  the  PHD  program,  may  be  manipulated  and 
changed.  If  a  new  DATASET  is  created  with  the  same  name  as  one  already  on 
mass  storage,  it  is  stored  as  a  new  version  of  the  file  in  systems  where  the 
operating  system  allows  this  (caution  must  be  exercised  when  using  multiple 
versions  as  a  purge  of  the  data  files  may  delete  all  but  the  last  version  of  a 
DATASET).  In  operating  systems  where  multiple  versions  are  not  allowed,  it 
replaces  the  old  version. 
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Most  of  the  aodulea  of  PHD  are  used  to  define  the  dataset  structure  (the 
base),  or  aanipulate  the  data  in  it.  In  figure  1,  on  the  following  page,  is  a 
typical  DATASET.  The  underlined  text  is  not  part  of  the  dataset  but 
explanatory  material. 


DATASET  NAME:  PRESUR 

(The  file  name  would  be  PRESUR. DAT ;n — n  is  the  version) 

COMMENT(one-line) :  7- 

-JAN-8 1/14: 02 

:03l VAPOR  PRESSURE 

FOR  PERWAXY  CHLORIDE 

Data  matrix: 

PSN:ID 

Pressure 

Temperature 

♦  Variable  NAME 

[Torr] 

[Kelvin] 

^  UNITS  of  variable 

0.1 

.001 

METRIC  of  variable 

2:Pnt  1 

34.52 

279.333 

♦  Observations  follow 

2:Pnt  2 

• 

16.6 

265.24 

• 

• 

• 

l:Pnt  99 

295.33 

• 

456.899 

• 

s 

*■  Last  observation 

Figure  1.  Sample  Data  Set 

In  the  example  DATASET  given  above  there  are  two  variables  stored  for  each 
observation.  The  operation  of  PHD  will  be  explained  in  terms  of  entering  the 
DATASET  given  above  and  performing  various  manipulations  on  it,  A  summary  of 
the  various  modules  and  the  information  (arguments)  they  may  require  or 
optionally  take  is  given  at  the  end  of  this  report. 

Hardware  requirements 

A  computer  capable  of  running  DEC's  BASlC-11  under  RSX-llM  or  VAX/VMS 
native  mode  BASIC  under  VMS  is  necessary  to  run  PHD  without  modification.  A 
mass  storage  device  capable  of  storing  the  DATASET  produced  and  sufficient 
main  siesiory  or  virtual  array  capability  as  defined  for  BASIC-11  is  also 
required.  The  Tektronix  A0I4  terminal  must  be  used  to  take  advantage  of  all 
the  features  of  PHD.  If  a  different  terminal  is  used,  spurious  output  of  8's, 
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9's,  r's  or  ;'8  will  occur  whenever  s  character  size  change  is  requested  by 
PBD.  The  .PLT  files  generated  by  PHD  under  RSX-llN  or  the  plots  produced  on 
the  VAX  nay  only  be  plotted  on  the  Tektronix  401A  or  a  .drminal  compatible 
with  it. 
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OPERATION 

Initiating  PHD 

Calling  the  PHD  program  varies  vith  operating  system  of  the  host 
computer.  Under  DEC's  RSXII-M  or  on  the  MINC,  PHD  is  called  as  any  other 
BASIC  program.  On  the  VAX  it  is  most  efficiently  called  as  a  program  under 
the  DCL  monitor.  Familiarity  with  the  computer  system  sufficient  to  sign  on 
and,  if  necessary,  get  into  BASIC  is  assumed. 

On  the  PDP  11/10  or  11/45  a  one  statement  BASIC  program  written  exactly  as 
follows  is  necessary  to  initiate  PHD.  Assuming  the  programs  for  PHD  arc 
stored  and  accessible  under  the  [142,1]  user  code  the  program  would  be  entered 
as  follows: 

NEW  PHm  (See  note  below) 

10  CHAIN  "[142,1] PHD"  LINE  lOO 

SAVED 

Once  such  a  program  is  written  and  saved,  PHD  may  be  initiated  at  any  time 
while  in  BASIC  by  entering  the  following: 

RUN  PHDD 

On  the  VAX,  the  PHD  program  is  used  in  a  compiled  form  normally.  The 
program  used  in  the  Chemistry  Directorate  (NC)  is  stored  in  the  directory 
[SEILERV.PHD]  as  PHD. EXE.  It  can  be  run  from  any  other  directory  by  using  the 
full  file  specification  in  the  command: 

RUN  [SEILERV.PHD] PHD 

Note:  Here  and  in  the  rest  of  this  report,  the  symbol  □  will  be  used  to 
represent  a  carriage  return. 
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By  using  the  DCL  global  definition, 

PHD:»RUN  (SEILERV.PHD]PHD 

the  program  may  be  initiated  simply  by  typing  in  PEDC  (in  all  NC  directories 
the  global  definition  is  made  in  the  login  comnand  file). 

From  the  initiation  of  the  PHD  program  until  the  user  exits  from  it,  the 
operations  appear  the  same  to  the  user  on  any  system.  To  continue  the 
illustration  of  using  PHD,  let  us  assume  the  user  proceeds  by  one  of  the  above 
methods  to  initiate  PHD. 

The  response  from  PHD  will  be  to  go  to  a  new  page  (clear  the  screen)  and 
present  the  MENU.  When  the  MENU  is  presented,  at  least  two  letters  of  one  of 
the  options  must  be  typed  in,  followed  by  any  argument  the  option  requires  or 
will  accept  (at  least  one  blank  must  separate  the  option  from  the  argument). 
For  example,  if  the  dataset  of  figure  1  is  to  be  created,  the  option  specified 
would  be  BASE,  and  the  argument  to  indicate  a  new  dataset  is  NEW.  The  comnand 
line  would  be  as  follows: 
base  NEWD 

PHD  responds  by  deleting  the  current  WORKSET  and  its  BASE  definition.  It 
then  requests  the  NAME  for  the  new  WORKSET  and  a  new  COMMENT  line.  New 
VARIABLE  definitions  may  then  be  inserted  into  the  BASE,  which  is  initially 
null  (i.e.,  no  variables  are  defined).  A  VARIABLE  is  added  by  specifying  its 
NAME, UNIT, METRIC  with  commas  between.  A  null  string  may  be  used  for  any  of 
the  three  fields  in  the  entry.  The  null  string  is  interpreted  to  mean  do  not 
change  from  the  current  value.  For  new  VARIABLES  the  initial  values  are  the 
null  string  for  NAME  and  UNIT  and  a  zero  for  the  METRIC.  A  zero  METRIC  means 
the  measurements  are  to  be  considered  exact.  For  the  example  DATASET  the 
following  entry  would  be  typical:  (Computer  responses  are  underlined.) 
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W0RKSET7  PRESURD 

COHMENT?  VAPOR  PRESSURE  FOR  PERWAXY  CHLORIDEO 
THE  CURRENT  BASE  HAS  0  VARIABLES  PER  OBSERVATION. 

MAME.UHIT, METRIC  FOR  ADDITIOHAL  VARIABLE  ?  PRESSURE, TORR, . In 
THE  CURREMT  BASE  HAS  1  VARIABLES  PER  OBSERVATION. 

PRESSURE  [TORR]  •»/-.! 

NAME, UNIT, METRIC  FOR  ADDITIONAL  VARIABLE  ?  TEMPERATURE, K, .00 ID 

THE  CURRENT  BASE  HAS  2  VARIABLES  PER  OrSERVATION. 

PRESSURE [TORR ]»/-.!  TEMPERATURE [k]*/-. 001 

NAME, UNIT, METRIC  FOR  ADDITIONAL  VARIABLE  ?  □ 

(/hen  the  BASE  has  been  defined,  one  may  enter  the  observations  using  the 
option,  ADD.  The  following  entry  provides  an  example: 

ADDO 

ENTERING  'ERROR'  WILL  CAUSE  REENTRY  OF  LAST  OBSERVATION  STARTED. 

NEXT  OBSERVATION  (/ILL  BE  NUMBER  1 

OBSERVATION  ID  OR  'END* ?2;PNTln 
PRESSURE [TOWl]?  34.1^ 

TEMPERATURE [K]?~ 279.3330 

NEXT  OBSERVATION  V/ILL  BE  NUMBER  2 

OBSERVATION  ID  OR  'END'? 

Succeeding  observations  would  be  entered  in  the  same  manner.  When  all 
observations  have  been  entered,  the  response  ENOn  is  given  in  response  to 
the  next  observation  ID  request. 

NEXT  OBSERVATION  WILL  BE  NUMBER  n 
OBSERVATION  ID  OR  *EWD'?  ENDO 

This  entry  will  terminate  the  ADD  module  and  return  control  to  the  MENU. 
Usually  at  this  point  it  is  prudent  to  safeguard  the  data  by  SAVEing  it,  or 
creating  a  DATASET  on  the  disk.  The  response  to  the  MENU  would  be  as  follova: 
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OPTIOH  T  SAVED 

CDRREHT  WOBKSET  WAKE:  <nork5et  nai^> 

HEW  DATASET  NAME  (cr  TO  USE  WORKSET  WAKE)?  <datasetaaiBe> 

The  WORKSET  would  then  be  stored  on  mass  storage  under  the  <datasetname>. 

As  can  be  seen  from  the  above  detailed  examples,  the  operation  of  the 
various  modules  involves  a  conversational  interaction  with  the  user.  For 
example,  the  response  required  b>  the  SAVE  module  is  for  a  DATASET  name  to 
store  the  data  under.  As  the  program  takes  some  action,  it  responds  with 
infonsational  replies  as  to  what  it  is  doing.  When  the  DATASET  has  been 
created  on  the  disk,  control  returns  to  the  MENU. 

After  entering  and  storing  the  data  in  a  DATASET,  a  listing  could  be  made 
to  check  for  errors.  The  option  LIST  will  present  the  DATASET  name,  COMMENT, 
and  data  matrix  listing  on  the  terminal  much  as  it  is  given  in  figure  1.  If 
errors  are  present,  the  CORRECT  module  may  be  used  to  make  changes.  The  data 
may  then  be  sorted  using  the  SORT  option,  PRINTed  on  the  line  printer,  FlTted, 
PLOTted  or  otherwise  manipulated  using  the  USE  option  to  call  programs  which 
have  been  written  by  the  user  in  BASIC. 

On  the  following  pages  are  details  on  the  use  of  each  of  the  options 
available.  Two  lists  of  the  BASIC  variables  used  internally  in  PHD  follow, 
one  for  the  RSX-llM  version,  the  other  for  the  user  accessible  BASIC  variables 
in  the  VAX/VMS  version.  Normally  the  BASIC  variables  used  by  PHD  should  be 
treated  as  read  only  variables.  If  it  is  necessary  to  change  the  values  of 
some  of  these  variables,  it  should  be  done  with  extreme  care  to  avoid 
disastrous  results  in  processing  data  in  a  user  module. 

TVo  short  descriptions  of  the  files  associated  with  particularly  useful 
PHD  options,  USE  and  FIT  are  then  given.  The  first  treats  the  writing  of  a 
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user  option  program  siodule  and  the  second  the  parasieter  description  files  used 
in  fitting  data  to  a  generalized  polynomial. 

Finally,  listings  of  PHD,  the  USER  envelope  program  ani  an  example  of  Its 
use  are  included.  A  summary  of  the  PHD  options,  with  their  arguments 
completes  the  report. 
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Detailed  Operation  Descriptions 
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Option;  ADD 
Syntax;  ADO 

nfaere  no  arguments  are  allowed. 

Description  of  results: 

Upon  entering  this  module  a  check  is  made  to  see  if  a  WORKSET  BASE  has 
been  defined  (i.e.t  if  the  number  of  variables  defined,  NZ,  is  greater  than 
zero).  If  it  has  not,  the  user  is  advised  to  use  the  modules  BASE  or  GET  to 
define  the  BASE  and  an  exit  is  made  to  the  MENU.  If  a  data  base  has  been 
defined,  the  following  advisory  is  given: 

ENTERING  'ERROR*  WILL  CAUSE  REENTRY  OF  LAST  OBSERVATION  STARTED. 

NEXT  OBSERVATION  WILL  BE  NUMBER  n 

OBSERVATION  ID  OR  'END*  ? 


One  of  the  following  valid  entries  is  made  at  this  point: 

□  which  is  the  same  as  entering  1 :□ 

n:D  where  n  is  a  valid  PSN  (plot  symbol  number),  se>.  PLOT 

option  for  valid  numbers.  No  check  is  made  of  the 
validity  at  this  point, 

:  S3  where  s  is  any  string  of  six  or  fewer  characters 

(comma  is  not  an  allowed  character), 
n:8a  where  s  and  n  are  as  defined  above  and  the  additional 

restriction  that  the  total  number  of  characters  in  the 
identifier  cannot  exceed  eight,  including  the  colon, 
ERRORP  this  entry  will  cause  the  previous  observation  to  be 
deleted  and  re-entry  of  that  observation  begun. 

ENDO  this  entry  will  terminate  the  ADD  option  and  return  to 

the  MENU. 

After  the  identifier  has  been  entered,  the  module  will  request  each  of  the 
variables  defined  in  the  data  base,  by  NAME,  in  the  appropriate  UNIT,  as 
follows: 

<nS3ie>[<unit>  ]  ? 


Only  three  possible  entries  are  valid  here: 

ENDO  this  entry  will  cause  an  EXIT  from  the  ADD  module  and 

the  partially  entered  data  observation  will  not  be 
included  in  the  dataset. 

ERROR3  this  entry  will  cause  the  current  observation  being 
entered  to  be  discarded  and  re-entry  begun. 
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Option;  ADD  (cont'd) 

<valueX3  vhere  <value>  is  a  real  nunber 

acceptable  to  BASIC,  a  positive  or  negative 
ninaber  in  the  range  10'**  to  10’*.  Any 
nunber  of  significant  digits  nay  be  entered, 
but  only  6  are  retained  accurately. 

Each  variable  will  be  requested  in  the  same  manner  until  all  are 
obtained.  No  more  than  1000  items  (number  of  variables  times  nuid>er  of 
observations)  may  be  entered  without  modification  of  PHD. 

Possible  errors  and  recovery  options: 

An  erroneous  value  may  be  re-entered  by  using  the  ERROR  entry  or  by 
changing  it  later  with  the  CORRECT  option.  Any  non-maeric  entry  for  a 
^alae>,  other  than  ERROR  or  END,  may  in  sobk  versions  result  in  an 
abnormal  exit  from  PHD  and  the  loss  of  the  current  NORXSET. 

Date  of  last  module  update:  23  Mar  82 
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Option;  BASE 

Syntax  1)  BASE 

2)  BASE  NEW 

3)  BASE  <dataset> 

%ihere  <dataaet>  is  the  name  of  a  valid  PHD  dataset. 

Description  of  results; 

For  the  first  syntax,  the  module  BASE  will  be  called  and  the  current 
WORKSET  BASE  will  be  displayed  as  follows: 

THE  CURRpT  BASE  lUS  n  VARIABLES  PER  OBSERVATION. 

^ame  l>[<unit  !>]■«•/-  <metric  1>  <name  2>  [<unit~2>  ]-*-/-<metric  2> 

ADD  I DELETE  I  CHANGE  NAMeTuNITS [METRIC  OF  A  VARIABLE  T 

where  <naiiie  i>,  <unit  i>  and  <inetric  i>  are  the  NAME,  UNIT, 
and  METRIC  of  VARIABLE  i.  There  are  four  valid  entries  at  this  point: 

□  this  entry  will  cause  a  return  to  the  MENU  with  no 

further  changes  in  the  BASE. 

AD  this  entry  will  result  in  a  request  for  a  new  VARIABLE 

definition  of  NAME, UNIT, METRIC.  A  null  entry  is  valid 
for  any  of  these  fields  (successive  commas  indicate  a 
leading  null  string).  If  there  is  more  than  one 
observation  in  the  workset,  PHD  will  repack  the  data 
at  this  point  giving  the  advisory  message:  REPACKING 
DATA  (STRETCH). 

DO  this  entry  will  cause  the  request  of  the  NAME  of  the 

variable  to  delete  from  the  workset.  This  name  must 
be  in  the  current  BASE  definition.  If  there  is  more 
than  one  observation  in  the  workset,  PHD  will  repack 
the  data  at  this  point  giving  the  advisory  message: 
REPACKING  DATA  (SQUEEZE). 

CD  this  entry  will  locate  the  VARIABLE  definition  in  the 

current  BASE  and  result  in  a  request  for  a  new 
variable  definition  of  NAME, UNIT, METRIC.  A  null  entry 
is  valid  for  any  of  these  fields  (successive  commas 
indicate  a  leading  null  string),  and  will  leave  that 
value  unchanged. 

The  second  syntax  results  in  a  discarding  of  the  current  WORKSET 
(equivalent  to  a  DELETE  BASE  option)  and  a  request  for  new  VARIABLE 
definitions  of  NAME, UN IT, METRIC  (See  the  A  response  in  syntax  1).  Additional 
VARIABLES  nay  be  defined  up  to  sixteen.  A  null  response  will  cause  return  to 
the  MEND. 
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Option;  BASE  (cont'd) 

The  third  syntax  results  in  the  opening  of  the  DATASET  stored  under  the 
name  given  and  the  construction  of  a  BASE  identical  to  the  one  specified  in 
that  DATASET. 

Possible  errors  and  recovery  options: 

Most  errors  may  be  corrected  by  using  the  C  response  of  syntax  1  above. 
Date  of  last  module  update;  23  Mar  82 
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Option;  CM  (COMMENT  Modification) 

Syntax:  CM  <optional  editing  character><coiiiaent> 

where  <coment>  is  any  character  string  which  does  not  include  a 
carriage  return,  and  <optional  editing  character>  an ^  its  possible 
meanings  are  as  follows: 

1)  >  meaning  replace  the  rightmost  portion  of  the  current 
COMMENT  with  the  following  string.  The  replacement  is  made  on  a 
character  for  character  basis. 

2)  <  meaning  replace  the  leftmost  portion  of  the  current 
COMMENT  with  the  following  string. 

3)  -  meaning  add  the  following  string  preceding  the  current 
COMMENT. 

4)  +  meaning  add  the  following  string  to  the  end  of  the  current 
COMMENT. 

Description  of  results: 

The  first  character  of  the  argument  is  examined,  if  it  is  an  editing 
character  the  current  COMMENT  is  modified  as  specified  by  that  editing 
character.  If  no  editing  character  is  found,  the  <comment>  becomes  the 
new  COMMENT  and  the  old  COMMENT  is  discarded. 

Possible  errors  and  recovery  options:  None 

Date  of  last  module  update:  27  Dec  82 
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Option:  CORRECT 

Syntax:  CORRECT 

where  no  arguments  are  allowed. 

Description  of  results: 

PHD  will  ask  for  the  sequence  number  of  the  observation  to  be  corrected: 
SEQUENCE  NUMBER  OF  OBSERVATION  TO  CORRECT  ?  nP 

where  n  is  the  sequence  numbe.  (obtained  from  the  LIST  or  PRINT  options). 
If  a  null  entry  (carriage  return  only)  is  given  here,  an  exit  is  made  to  the 
main  MENU. 

Having  received  the  sequence  number  of  an  observation,  PHD  checks  it  to  make 
sure  it  is  in  range  (see  possible  errors).  The  observation  is  then  presented 
on  the  terminal  and  a  request  for  corrections  made. 

'PSN:ID*  OR  ’VARIABLE  NAME, VALUE*  ? 

There  are  several  possible  entries  at  this  point: 

Return  to  get  new  sequence  number  of 
observation  to  be  corrected. 

This  will  result  in  only  the  plot  symbol 
number  being  changed. 

This  results  in  only  the  identifier  being 
changed. 

This  changes  both  the  plot  symbol  number  and 
the  identifier.  Any  of  the  options 
containing  a  colon  (:)  will  result  in  the 
message  NEW  ID  IS  <psn>:<id>  and  a 
request  for  further  corrections  for  the  same 
observat ion . 

<variable  natne>,<valuexi 

This  results  in  the  value  of  the  variable 
named  being  changed  to  <value>  for  the 
selected  observation.  After  the  change  has 
been  made,  requests  for  further  changes  to 
this  observation  will  be  made. 

Possible  errors  and  recovery  options: 

The  presence  of  the  colon  (:)  in  the  input  indicates  that  entry  is  a 
change  to  either  the  plot  symbol  number  or  identifier,  depending  on  the 
location  of  the  colon.  No  check  is  made  on  the  validity  of  the  plot  symbol 
number  or  the  identifier.  The  new  plot  symbol  number  and  identifier  are 
displayed  after  the  change.  Note  that  the  composite  string  <p8n>:<id> 
is  truncated  after  eight  characters  with  no  warning. 


a 

<psn> :□ 

;<id>0 

<psn>:<id>D 
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Option:  CORRECT  (Cont'd) 

If  the  input  does  not  correspond  to  any  of  the  above  syntaxes,  then 
UHACCEPTABLE  INPUT  is  printed  and  the  program  returns  to  the  request  for  a 
sequence  number. 

In  the  last  syntax,  if  a  variable  name  is  referenced  that  does  not  exist 
in  the  current  BASE,  a  warning  is  printed,  NO  VARIABLE  <variable  name> 

FOUiro,  no  chafes  are  made  to  the  observation,  and  the  next  correction  for 
that  observation  is  requested. 
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Option:  DELETE 

Syntax:  1 )  DELETE 

2)  DELETE  <sequence  8et> 

where  <sequence  8et>  is  a  set  of  <sequence  range>s  separated 
by  commas  ( , ) . 

<sequence  range>  is  <sequence  nuinber>  or  <start 
sequence  number>-  or  <start  sequence  nuiiiber>-<end 
sequence  nunber>  or  — <end  sequence  nuaber> 

Description  of  results: 

If  no  <8equence  set>  is  specified,  as  in  syntax  1,  then  PHD  will 
request  one  with  Che  message: 

SEQUEHCE  NUMBER(S)  OF  OBSERVATION(S )  TO  DELETE  ? 

PHD  will  scan  the  <sequence  set>  entered  and  print  the  message: 

MARKED  FOR  DELETION 

ID  variablename  l[unitl  variablename  2 [unit]  etc . 

(A  listing  of  the  observations  marked  will  follow  this  header  message.) 

The  request  for  a  <sequence  set>  given  above,  will  also  occur  each 
time  after  all  observations  in  a  <sequence  set>  have  been  marked  for 
deletion.  Another  sequence  set  may  then  be  entered  or  a  null  entry  may  be 
made  in  order  to  return  to  the  MENU.  When  a  return  to  the  menu  is  made,  PHD 
deletes  all  the  observations  that  have  been  marked  for  deletion  and 
resequences  the  observations  that  are  left,  issuing  the  message 
n  OBSERVATION(S)  DELETED;  WORKSET  RESEQUENCED,  on  exit  to  the  MENU. 

Possible  errors  and  recovery  options: 

1)  If  a  <sequence  range>  is  illegal,  (i.e.,  one  in  which  <8tarting 
sequence  number>  is  larger  chan  <end  sequence  number>)  then  the  message 
ILLEGAL  SEQUENCE  SET  is  printed,  no  observations  are  marked  for  deletion  and 
the  scan  of  the  <sequence  set>  continues. 

2)  If  a  <8equence  range>  extends  beyond  the  current  set  of  data 
sequence  number,  i.e.,  less  than  1  or  greater  than  the  number  of  observations, 
then  the  message  SEQUENCE  RANGE  IS  FROM  1  TO  n  where  n  is  the  sequence  number 
of  the  last  observation. 

3)  If  a  observation  already  marked  for  deletion  occurs  in  an  another 
<sequence  range>  then  the  message  OBSERVATION  #  n  ALREADY  MARKED  FOR 
DELETION  will  be  issued.  The  observation  will  remain  marked  for  deletion. 

Date  of  last  module  update;  27  Dec  82 
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Option:  ENTER 
Syntax;  ENTER 

where  no  arguments  are  allowed. 

Description  of  results: 

The  current  WORKSET  BASE  description  will  be  printed  in  the  format  of  the 
following  example: 

NAME  UNIT  METRIC 

ID 

TEMP  DEG  C  0 

SP  POND  MHO/CM  .0002 


SEQ  NR 
0 
1_ 

2 


A  request  for  which  VARIABLE  is  to  be  entered  into  the  WORKSET  will  be 

made.  In  the  following  example  the  VARIABLE  "SP  COND"  is  to  be  entered. 

NUMBER  OF  VARIABLE  TO  ENTER  ?  2D 

An  advisory  message  and  request  for  change  for  the  specified  variable  in 
each  observation  will  be  made  as  in  the  following  example: 

ENTER  NEW  VALUES  FOR  SP  COND 

i  2:pt  1  .59687  NEW  VALUE  (CR  FOR  NO  CHANGE)>  ?  .5968to 

I  2:pt  2  .A7859  NEW  VALUE  (CR  FOR  NO  CHANGeT>  ?  n 


After  the  last  value  of  the  VARIABLE  is  entered,  a  return  will  be  made  to 
the  main  MENU. 

A  different  format  is  used  if  the  ID  is  to  be  changed.  The  initial 
display  after  entering  the  value  0  in  response  to  the  number  of  the  variable 
to  enter  is  as  follows: 

ENTER  NEW  ID'S 

IF  ALL  NEW  ID'S  ARE  TO  BE  ALIKE,  ENTER  NEW  ID  HERE  (OTHERWISE  CR)>  ? 

If  any  of  the  ID's  are  to  be  different,  then  all  of  them  must  be  entered 
individually  (or  left  as  they  are),  unless  the  user  wishes  to  change  them  all 
and  then  make  selective  changes  by  using  ENTER  again.  If  a  non-null  entry  is 
made  at  this  point,  then  all  ID's  are  changed  to  that  value  and  the  advisory 
message  DONE  is  printed  and  an  exit  is  made  to  the  main  MENU.  If  a  null  entry 
is  made  here  then  the  changes  made  are  much  the  same  as  for  numerical  values. 
The  display  of  the  individual  observation  ID's  for  change  appears  as: 


21 


PHD:  Version  2.2 


Option:  ENTER  (Coat ’d) 

OBSERVATION  ni  PRESENT  PSH:ID  IS  ENTER  HEW  PSW:ID  (CR  FDR  NO 

CHANGE)?  8:NEWIin 

tfhere  n  is  the  observation  sequence  number  and  mriiiiii  is  the  current 
<psn>:<id>  for  observation  n.  In  the  example,  the  new  PSN  and  ID 
would  be  8  and  NEUIO,  respectively. 

Possible  errors  and  recovery  options: 

Please  note  in  changing  ID'r,  that  both  the  PSN  and  the  ID  must  be  entered 
separated  by  a  colon  (:).  I£  the  changed  PSN: ID  string  does  not  have  a  number 
in  the  range  1-12  followed  by  a  colon  as  its  initial  characters,  erroneous 
results  may  be  obtained  in  SAVEing  the  data  and  in  the  PLOT  module  as  well  as 
other  modules  which  depend  on  this  specific  format  for  the  PSN:ID  string. 

Note  also  that  the  total  length  of  the  PSN: ID  string  including  the  colon  must 
be  less  than  9  characters.  If  the  string  is  changed  to  one  that  is  longer,  it 
will  be  truncated  after  the  first  eight  characters  with  no  warning  message. 

If  a  VARIABLE  sequence  number  that  is  out  of  the  range  for  the  current 
RASE,  no  warning  message  is  given;  instead,  the  BASE  description  is 
redisplayed  and  another  request  for  the  number  of  the  VARIABLE  to  enter  will 
be  made. 

If  an  attempt  is  made  to  replace  an  observation  value  with  a  character 
string  that  is  not  a  valid  BASIC  number,  ILLEGAL  NUMBER  will  be  displayed,  and 
a  new  request  for  a  change  in  that  observation  made. 

Date  of  last  module  update:  27  Dec  82 
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Option:  EXIT 

Syntax:  DEXIT 

2)  EXIT  < program  file  naine> 

where  <program  file  namE:>  ia  the  name  of  a  valid  .EXE  file  to  which 
control  ia  to  be  transferred  on  exit  from  PHD. 

Description  of  results: 

PHD  will  exit  to  the  program  file  if  specified  (syntax  2)  or  back  to  the 
operating  system  or  BASIC,  whichever  it  was  called  from  (syntax  1).  Note 
that  on  leaving  PHD  all  data  currently  in  the  UORKSBT  is  lost.  If  the  data 
has  not  previously  been  SAVEd  in  a  peraanent  DATASET,  it  should  be  saved 
before  exiting  by  using  the  SAVE  option. 

Possible  errors  and  recovery  options:  None . 

Date  of  last  module  update:  27  Dec  82 


I 

i 
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Option:  FIT 

Syntax:  1)  FIT 

2)  FIT  <paraiDeter  file  name> 

where  <parameter  file  name>  is  a  name  for  a  file  that  provides 
equation  and  parameter  information  to  fit  to  the  observations  in  the  current 
WORKSET.  The  file  extension  must  be  .PAR  (see  Parameter  Files  for  Data 
Fitting,  p  38) 

Description  of  results: 

If  syntax  1  is  used,  PHD  will  ask  for  the  parameter  file  name  by 
requesting: 

NAME  OF  PARAMETER  FILE  ?  . 

After  the  parameter  file  is  entered  results  are  the  same  as  for  syntax  2. 

With  syntax  2,  PHD  will  read  the  parameter  file  specified.  The  banner 
(fit  title)  is  printed  as  it  is  read  and  an  indication  of  end-of-file  is  given 
when  the  end  of  the  file  is  reached.  The  parameters  specified  are  then  sorted 
based  on  their  DECODE  symbol.  D‘. agnostics  are  printed  as  the  parameters  are 
checked.  A  typical  output  stream  follows: 

OPTION?  FIT  IRATE 

LINEAR  FIT  OF  %  YIELD  TO  TIME  (INIT.  RATE) 

**]e6f** - 

3  PARAMETERS.  2  ADJUSTABLE. 

SORTING. . . 

PARAMETER  TABLE 


NR  NAME 

UNITS 

EQ  DECODE 

INIT  VALUE 

ADJ 

1  INT 

IHHH 

<+00> 

0 

Y 

2  IRATE 

(ZPROD/MIN) 

<  +  10> 

0 

Y 

3  IPROD 

<-01> 

-1 

N 

After  the  parameters  are  checked,  two  inputs  are  requested:  A  new  banner 
and  a  criterion  for  fit  convergence.  Default  values  appear  in  brackets  in  the 
prompts . 

BEGIN  FIT: 

TITLE:  [LINEAR  FIT  OF  %  YIELD  TO  TIME  (INIT.  RATE)]?  %YIELD  FOR  ACETONED 

DELTA  SIGMA  FOR  CONVERGENCE  [1E-7|?0 

The  current  banner  is  displayed  (taken  from  the  parameter  file)  and  if  a 
different  one  is  desired,  it  should  be  entered  after  the  question  mark.  A 
null  entry  will  retain  the  current  banner.  The  delta  sigma  requested  is  the 
upper  bound  for  the  fractional  change  in  standard  deviation  which  will 
terminate  the  fitting  process.  Any  fractional  change  in  the  standard 
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Option:  PIT  (Cont'd) 

deviation  less  than  the  specified  delta  sigma  will  satisfy  the  criterion  for 
convergence.  During  the  fitting  process  changes  in  the  adjustable  paraaieters , 
sigma  and  delta  sigma  are  printed  after  each  iteration. 

ITERATION  1  S1Q<A  =  131.789 
7  OBSERVATIONS  USED.  DELTA  SIGMA  =  0 

ADJUSTABLE  PARAMETERS  CORRECTED  VALUE 
INT  =  132.385  IRATE  -2.31018 

ITERATION  2  SICMA  -  .125^67 
7  OBSERVATIONS  USED.  DELTA  SIGMA  =  .983 

ADJUSTABLE  PARAMETERS  CORRECTED  VALUE 
INT  «=  132.385  IRATE  -2.31 

ITERATION  2  CONVERGENCE  OCCURRED. 

HIT  CARRIAGE  RETURN  TO  CONTINUE.  ? 


After  the  indication  chat  the  fitting  process  has  converged,  a  pause  is 
made  to  await  a  carriage  return  from  the  user.  On  recieving  this  carriage 
return,  the  screen  will  be  cleared  and  the  name  of  the  dataset  fitted,  the 
banner,  sigma  (the  standaro  deviation)  and  delta  sigma  and  a  sunmiary  of  the 
parameters  and  their  uncertainty  is  given.  Continuing  with  the  above  example: 

CPY26BF 

XYIELD  FOR  ACETONE 


SIGMA  =  .125667  DELTA  SIGMA  =  -1.99E-8 


PARAMETER 

UNITS 

VALUE  AND  UNCERTAINTY 

INT 

l%] 

(+/-)  8.5E-02 

IRATE 

[ZPROD/MIN] 

(+/-)  3.2E-01 

fPROD 

rz] 

NOT  ADJUSTABLE 

FIT  TERMINATION-ENTER  cr  FOR  MAIN  MENU 
- 


On  display  of  the  FIT  TERMINATION  message  the  program  will  wait  for  a 
carriage  return  before  exiting  to  the  main  MENU. 

Possible  errors  and  recovery  options : 

If  the  specified  parameter  file  does  not  exist  the  message  "NO  SUCH  FILE” 
will  be  given  and  a  the  module  reinitiated. 

There  is  a  maximum  of  16  parameters  (adjustable  and  non-ad justable 
combined).  If  an  attempt  is  made  to  enter  more  than  that,  the  message  MAXIMUM 
P ARAMETERS  EXCEEDED  will  be  printed  and  an  attempt  to  fit  will  be  made  with 
the  parameters  already  specified. 
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Opt ion;  FIT  (Cont'd) 

During  the  parameter  checking  phase  a  number  of  different  messages  siay  be 
given: 

If  no  decode  string  is  given  for  a  parameter,  HAS  HULL  DECODE  is  printed 
and  a  DECODE  of  <-0>  is  assumed  (ASSUMING  DECODE  <-0>) . 

If  the  DECODE  does  not  specify  the  exact  number  of  VARIABLES  in  the 
current  WORKSET,  then  one  of  two  messages  is  given: 

1)  if  it  specifies  fewer  than  the  number  of  VARIABLES  in  the 
WORKSET  then  ACCESSES  VARIABLES  n-m  is  printed. 

2)  if  it  specifies  more  than  the  number  of  VARIABLES  in  the 
WORKSET,  then  DECODE  SPECIFIES  NON-EXISTENT  VARIABLES  j-k 

DECODE  CHANGED  FROM  <current  decode>  TO 
<truncated  decode>. 

If  an  illegal  character  is  present  in  the  DECODE,  it  is  replaced  by  a  0 
and  the  messages  USES  ILLEGAL  DECODE  ELEMENT  '< illegal  char>' 
and  ILLEGAL  ELEMENT  REPLACED  BY  *0',  are  given. 

Date  of  last  module  update:  27  Dec  82 
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Option:  GET 

Syntax:  1)GET 

2)GET  <data8et> 

where  <data8et>  is  the  name  of  an  existing  DATASET  on  the  mass  storage 
device  (the  file  name  is  <dataset>.DAT) ,  The  format  of  the  DATASET  file 
on  mass  storage  must  be  compatible  with  PHD. 

Description  of  results: 

If  syntax  1  is  used,  PHD  will  repeatedly  request  <data8et>s  and  add 
the  observations  from  the  specified  DATASET  to  the  current  WORKSET  until  a 
null  entry  is  made  for  the  <data8et>.  The  BASE  definition  of  the  first 
DATASET  placed  in  the  WORKSET  will  be  retained  when  additional  DATASETs  are 
added  to  the  WORKSET.  When  a  null  entry  is  given  an  exit  to  the  main  MEMO  is 
made.  After  the  initial  DATASET  is  read  in  ,  each  DATASET  is  checked  for  the 
number  of  VARIABLES  per  observation,  (N2).  If  the  number  of  VARIABLES  per 
observation  is  not  the  same  as  in  the  WORKSET,  the  DATASET  is  not  added  to  the 
WORKSET  and  an  advisory  message  is  printed  (MISMATCH  IN  NUMBER  OF  VARIABLES 
PER  OBSERVATION).  The  COMMENT  and  variable  NAMEs ,  UNITS,  and  METRICS  are 
ignored  for  any  DATASET  if  a  WORKSET  BASE  already  exists  and  the  original 
value  of  these  items  is  retained. 

When  syntax  2  is  used,  only  the  DATASET  specified  is  added  to  the  WORKSET, 
then  an  exit  is  made  back  to  the  MENU. 

In  both  syntaxes,  a  number  of  advisory  messages  are  printed  while  the 
DATASET  is  being  processed: 

DATASET:  <dataset  name>  n  VARIABLES  PER  OBSERVATION. 

TOMMENT  line 

(a  summary  of  the  BASE  definition  is  printed  here) 

On  completion  of  the  DATASET  processing,  the  following  messages  are 
printed : 

NOW  nn  OBSERVATION( S )  IN  CURRENT  WORKSET 

DATASET  <data8et  name>  CLOSED. 


Possible  errors  and  recovery  options: 

If  a  <data8et>  that  does  not  exist  on  mass  storage  is  specified,  the 
message  NO  SUCH  DATASET  is  given.  Badly  formatted  data  may  produce  the  error 
messages  WO  DATA  IN  OBSERVATION>  pen: id .  or  ERROR  nn  IN  LINE  mmm  IN  GET. 

Tn  the  former  case  undefined  (and  probably  erroneous)  data  will  exist  in  the 
observation  specified.  The  second  error  is  more  serious  and  will  terminate 
the  GET  operation  at  the  point  where  the  error  occurred  with  a  return  to  the 
main  MENU.  In  this  case  the  state  of  the  WORKSET  will  be  unknown. 

Date  of  last  module  update:  27  Dec  82 
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Option:  LIST 
Syntax:  LIST 

where  no  arguments  are  allowed. 

Description  of  results: 

The  contents  of  the  current  WORKSET  are  listed  on  the  terminal  in  a 
tabular  form  along  with  the  WORKSET  name  and  CCMMENT.  If  more  than  30  items 
are  in  the  WORKSET,  the  listing  is  halted  every  50  entries  for  viewing, 
copying,  etc.  Entering  a  carriage  return  will  resume  the  listing  process. 

Possible  errors  and  recovery  options:  None . 

Date  of  last  module  update;  27  Dec  82 
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Option:  MFP  (Indirect  Option  File  Processing) 

Syntax:  MFP 

where  no  arguments  are  allowed. 

Description  of  results: 

PHD  will  request  the  name  of  a  command  file,  <@file>.  This  command 
file  must  contain  a  sequence  of  PHD  MENU  options.  The  file  name  on  mass 
storage  must  be  <@file>.PHD.  This  type  of  file  must  have  been  previously 
built  using  one  of  the  host  system  editors  (e.g.,  the  MCR  EDI  editor,  or  EDT 
under  DCL).  When  in  the  MFP  mode,  instead  of  returning  to  the  main  MENU  for 
instructions,  PHD  will  process  the  options  from  the  comnand  file  specified 
instead.  Note  that  sub-commands  (i.e.  those  given  in  response  to  requests 
within  each  module)  must  still  be  input  from  the  terminal  during  execution  of 
that  module.  The  occurence  of  an  EXIT  command  in  the  command  file  will  cause 
a  termination  of  the  MFP  mode  and  a  return  to  the  main  MENU  at  the  terminal. 

Possible  errors  and  recovery  options: 

Illegal  commands  and  arguments  are  treated  the  same  as  if  they  had  been 
submitted  in  response  to  the  main  MENU.  Failure  to  include  an  EXIT  at  the 
end  of  the  command  file  will  result  in  the  abnormal  termination  of  PHD  and 
loas  of  the  current  UORKSET. 
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Option;  PLOT 
Syntax:  PLOT 

where  no  arguments  are  allowed. 

Description  of  results: 

This  module  is  designed  to  plot  all  or  portions  of  the  current  WORKSET. 
Vfhen  called  the  PLOT  module  will  request  information  on  what  observations  to 
plot,  VARIABLES  and  labels  for  the  axes  and  information  necessary  for  "nice" 
scaling  and  plotting.  All  plots  are  made  in  the  same  format:  A  frame  with 
tic  marks  at  the  specified  intervals,  with  labels  on  the  left  of  the  ordinate 
and  below  the  abcissa.  A  title  is  placed  across  the  bottom  of  the  plot, 
centered  below  the  abcissa  label. 

The  plot  can  be  made  only  on  a  TEKTRONIX  40xx  series  terminal  or  one 
compatible  with  these  terminals.  An  example  of  the  dialog  occurring  in  a  call 
to  PLOT  follows.  In  the  example  below,  all  values  were  taken  as  the  default 
by  simply  entering  a  carriage  return.  If  desired,  changes  could  be  made  by 
entering  the  desired  value  (character  string  or  number).  : 

PLOTn 

WELCOME  TO  PhD  PLOT  II 
n  VARIABLES/OBSERVATION 
m  OBSERVATIONS 
VARIABLES  AVAILABLE: 

VARIABLE  #  NAME 


1  namel 

2  name2 

3  name3 


VARIABLE  AS  ABCISSA[1]  ?n 
VARIABLE  AS  ORDINATE [2]  ?n 
FIRST  OBSERVATION  TO  PL0T{11  ?□ 
LAST  OBSERVATION  TO  PI.OT[m]  ?n 
PLOTTING  namel  IN  unitl 
VERSUS  name2  IN  unit2 
OBSERVATIONS  1-m 


UNITS 

unitl 
unit2 
unit  3 


TO  LEAVE  UNCHANGED,  PRESS  <CR> 

PLOT  TITLE:  COMMENT  for  current  WORKSETO 


X 

LABEL; 

namel? 

UNITS: 

unitl? 

Y 

LABEL: 

namel? 

UNITS: 

unitl? 
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Option:  PLOT  (Cont'd) 

AXIS  LIMITS  (TO  LEAVE  UNCHANGED,  PRESS  <CR>) 


X  MINIMUM  ^  xmin?  O 
X  MAXIMUM  =  xmax?  Q 

Y  MINIMUM  =  ymin?  O 

Y  MAXIMUM  =  ymax?  O 

TIC  MARK  INTERVALS  (TO  LEAVE  UNCHANGED ,  PRESS  <CR>) 

X:  xtic?  n 

Y:  ytic?  □ 

MIN  X  =  xdatamin  MAX  X  =  xdatamax 

MIN  Y  =  -ydatainin  MAX  Y  =  ydatamax 

SYMBOL  SIZE  (1-10)  [Sj  ?U 

After  the  entry  for  the  symbol  size,  the  screen  will  be  cleared  and  the 
plot  produced.  PHD  then  awaits  a  carriage  return  to  again  blank  the  screen 
and  return  to  the  main  MEMtf. 

The  legal  plot  symbol  numbers  (PSN  in  the  PSN:1D  string  for  each 
observation)  and  the  symbols  they  produce  are  as  follows: 

1  + 

2  X 

3  * 

4  n 

5  cl  iamond 

6  & 

7  V 

8  Move  to  location  but  plot  no  symbol  (move  with  pen  up). 

9  Solid  line  from  current  location  to  specified  location 
(move  with  pen  down). 

10  Same  as  9  but  dotted  line. 

11  Same  as  9  but  dash-dot  line. 

12  Same  as  9  but  dashed  line. 

Possible  errors  and  recovoi ■  options: 

In  general,  the  PLOT  module  will  recover  with  an  error  message  ILLEGAL 
NUMBER  or  the  system  error  message  for  "number  required",  whenever  a 
non-numeric  is  entered  when  a  number  is  required.  The  labels  have  length 
limitations  which  depend  on  the  length  of  the  UNIT  name  and  scaling  factors 
for  the  particular  plot.  Generally  no  warning  message  is  gi'*en  when  a  label 
is  truncated,  except  its  truncated  form  on  the  plot  itsel'  . 

Date  of  last  module  update:  27  Dec  82 
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Option;  PRINT 

Syntax;  1) PRINT 

2)PRINT  <nuniber  of  copie8> 

where  <nuiiiber  of  copies>  is  an  integer  between  I  and  10. 

Description  of  results; 

A  listing  file  (PHD. LIS)  containing  the  number  of  copies  of  the  current 
WORKSET  specified  in  <nutnber  of  copies>  is  produced  on  the  mass  storage 
device.  The  file  may  be  spooled  to  the  printer  using  the  DCL  command  PRINT 
PHD  in  response  to  the  $  aftet  an  EXIT  has  been  made  from  PHD. 

Possible  errors  and  recovery  options:  None . 

Date  of  last  module  update:  27  Dec  82 
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Option;  RENAME 

Syntax:  1)  REJIAME 

2)  RENAME  <new  workset  naiiie> 

where  <new  workset  naiiie>  is  a  file  name  acceptable  to  VMS. 

Description  of  results: 

For  both  syntaxes,  PHD  will  display  the  current  WORKSET  name  inaaediately 
on  entry: 

CURRENT  WORKSET  NAME:  <work?et  name> 


If  syntax  1  has  been  used,  a  <new  workset  naine>  will  then  be  requested: 
WORKSET  NAME  ? 

The  WORKSET  name  will  be  changed  to  the  one  specified.  If  a  null  entry  is 
given,  the  WORKSET  name  will  not  be  changed  and  an  immediate  exit  to  the  main 
MENU  will  be  made. 

After  the  WORKSET  name  change  is  made,  the  current  COMMENT  line  will  be 
displayed: 

CURRENT  COMMENT:  <current  coniment> 


Changes  may  then  be  made  to  the  COMMENT  line  by  using  a  <coinnient> 
optionally  preceded  by  an  editing  symbol  (see  CM  option),  when  the  request, 

COMMENT  ? 

is  made.  If  a  null  entry  is  made  at  this  point,  the  COMMENT  will  remain 
unchanged.  A  return  is  made  to  the  main  MENU  after  the  request  for  COMMENT 
editing  has  been  processed. 

Possible  errors  and  recovery  options: 

A  WORKSET  name  longer  than  9  characters  may  cause  problems  in  SAVEing  the 
WORKSET.  COMMENTS  longer  than  80  characters  will  be  truncated  to  80 
characters  on  transfer  to  any  user  written  module  (USE  option).  COMMENT  lines 
longer  than  40  characters  may  be  truncated  in  the  title  of  any  plot  (PLOT 
opt  ion) . 
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Option;  SAVE 

Syntax:  1 )  SAVE 

2)  SAVE  <data8et  name> 
where  <data8et  naiiie>  18  a  valid  VMS  file  name. 

Description  of  results: 

If  syntax  1  is  used,  PHD  will  request  a  <dataset  na8ie>,  prompting  the 
current  HORKSET  naote  as  a  default. 

CURRENT  WDRKSET;  <workset  name> 

NEW  DATASET  NAME  (cr  IF  NO  CHANGE)  ?  <dataset  nanie>D 

A  null  return  at  this  point  will  result  in  the  WORKSET  being  written  as  a 
DATASET  with  the  same  name  as  the  current  WORKSET.  The  WORKSET  will  then  be 
written  to  mass  storage  under  the  <dataset  name>  specified  (if  a  DATASET 
with  that  name  already  exists  on  mass  storage  then  a  new  version  will  be 
created),  with  the  following  advisory  messages: 

CREATING  NEW  DATASET;  <dataset  name> 
n  VARIABLES  PER  OBSERVATION,  m  OBSERVATIONS. 

(a  display  of  the  BASE  being  written  out  will  appear  here) 

DATASET  <dataset  name>  CLOSED. 

An  exit  to  the  main  MENU  then  occurs. 

Possible  errors  and  recovery  options:  None . 

Date  of  last  module  update;  27  Dec  82 
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Option:  SORT 
Syntax:  SORT 

where  no  arguments  are  allowed. 

Description  of  results: 

The  sorting  routine  is  designed  to  implement  a  hierarchy  of  sorts.  The 
latest  sort  made  is  always  the  major  sort,  the  earliest  sort  the  most  minor. 
The  observations  will  be  in  the  order  specified  by  the  major  sort  VARIABLE 
first.  When  values  of  the  major  VARIABLE  are  identical,  the  order  will  be 
determined  by  the  value  of  the  VARIABLE  just  minor  to  the  major  VARIABLE,  etc. 
to  the  most  minor  sort  specified. 

On  entry  into  this  module,  a  sequenced  list  of  VARIABLES  in  the  WORKSET 
BASE  will  be  given  along  with  an  advisory  message: 

THERE  ARE  n  VARIABLES  PER  OBSERVATION. 

(sequenced  list  of  VARIABLES  appears  here) 

A  POSITIVE  SEQUENCE  NUMBER  IMPLIES  ASCENDING  SORT 

A  NEGATIVE  SEQUENCE  NUMBER  IMPLIES  DESCENDING  SORT 

'CR*  WILL  IMPLEMENT  THE  SORT  VECTORS  THAT  HAVE  BEEN  BUILT 

SEQUENCE  NUMBER  OF  VARIABLE  TO  SORT  ON? 

In  addition  to  the  sequence  number  for  a  VARIABLE,  ID  or  -ID  may  be 
entered  to  build  a  sort  vector,  alphabetically  or  in  reverse  alphabetical 
order,  respectively.  In  these  cases  the  sort  will  be  made  on  the  value  of  the 
<psn>:<id>  for  each  observation. 

As  each  VARIABLE  sequence  number  entry  is  made,  the  advisory  BUILDING  SORT 
VECTOR  ON  KEY  <variable  name> .  will  be  given  and  the  sort  vector  modified 
to  reflect  this  sort^  On  completion  of  the  vector,  VECTOR  BUILT  will  be  given 
and  a  request  for  another  VARIABLE  to  sort  on  made.  When  the  final  vector  has 
been  built,  the  actual  sort  (physical  rearrangement  of  data  in  memory)  is  made 
by  giving  a  null  entry  (carriage  return  only). 

The  sort  vector  built  to  that  point  is  then  implemented.  The  advisory 
SORT  IN  PROGRESS  is  given  and  a  period  (.)  is  printed  for  every  vector  segment 
processed  during  the  sort.  When  the  sort  is  done  the  message  SORT  COMPLETED. 
is  given  and  an  exit  to  the  main  MENU  made. 

Possible  errors  and  recovery  options: 

If  a  VARIABLE  sequence  number  that  is  out  of  range  is  given,  the  entry 
advisories  will  be  repeated  and  a  new  request  for  sequence  number  made.  An 
alphxnoaeric  entry  that  is  not  either  ID  or  —ID  or  nail  will  result  in 
nbnorwsl  termination  of  PHD  and  the  loss  of  the  current  flOBKSET. 

Date  of  last  module  update;  27  Dec  82 
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Option:  SWAP 
Syntax:  SWAP 

where  no  arguments  are  allowed. 

Description  of  results: 

This  option  is  useful  where  a  particular  ordering  of  the  VARIABLES  in  the 
WORKSET  is  required.  Options  such  as  the  FIT  option  or  user  written  options 
may  depend  on  having  the  VARIABLES  in  a  particular  order  within  an 
observation.  When  two  different  options  require  different  ordering  of 
VARIABLES,  it  becomes  necessary  to  either  create  two  DATASETs  or  reorder  the 
VARIABLES.  To  avoid  this  difficulty  the  SWAP  option  allows  the  exhange  of  two 
different  VARIABLES  within  all  observations.  VARIABLES  may  be  placed  in  any 
order  desired  by  a  suitable  sequence  of  SWAP  operations. 

On  entry  into  this  module,  a  sequenced  listing  of  the  VARIABLES  in  the 
BASE  will  be  given  and  a  request  for  the  sequence  numbers  of  the  VARIABLES  to 
be  swapped  made. 

VARIABLES  IN  FILE 

( sequenced  1 ist ing  of  variables  in  the  WORKSET  BASE  appears  here) 

NUMBERS  OF  VARIABLES  TO  EXCHANGE  ? 

Valid  sequence  numbers  for  two  VARIABLES  to  be  exchanged  must  be  given  at 
this  point.  The  values  of  the  VARIABLES  for  all  observations  in  the  WORKSET, 
their  NAMEs,  UNITs,  and  METRICS  will  be  exchanged.  A  return  to  the  main  MENU 
is  then  made. 

Possible  errors  and  recovery  options: 

If  invalid  numbers  are  entered  as  VARIABLE  sequence  numbers,  a  system 
error  will  result.  If  an  out-of-range  sequence  number  is  given,  no  action 
will  be  taken  and  the  SWAP  module  re-initiated. 

Date  of  last  module  update:  27  Dec  82 
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Option;  OSE 

Syntax:  1)  OSE 

2)  USE  <user  file  naiBe> 

where  <U8er  file  nane>  must  be  the  nane  of  an  .EXE  type  file  built 
using  the  USER.BAS  aodule  as  an  envelope  program  (See  How  to  Write  a  User 
Program  Module,  p  35).  The  user  written  program  lines  which  may  modify  the 
WORKSET  ouist  be  sequenced  between  1000  and  29998  in  order  not  to  overlap  with 
the  stateaients  of  the  envelope  program.. 

Description  of  results: 

If  syntax  1  is  used,  PHD  will  request  the  <user  file  name>: 

PROGRAM  NAME  ? 


If  a  null  entry  is  made  at  this  point  then  a  return  to  the  main  MENU  will 
be  made;  otherwise,  the  action  proceeds  as  for  syntax  2. 

When  syntax  2  is  used,  a  request  is  made  for  arguments  to  be  passed  to  the 
user  program  via  the  string  variable  Yj.  The  WORKSET  is  stored  temporarily  in 
the  file  USER.TMP  and  the  user  program  called.  When  the  user  program  returns 
to  PHD,  a  check  is  made  for  the  existence  of  the  file  USER.TMP.  If  this  file 
is  found  an  advisory,  USER  PROGRAM  COMPLETE. . .RECOVERING  WORKSET.  is  made  and 
the  contents  of  the  file  USER.TMP  are  read  into  the  WORKSET.  The  retrieval 
and  temporary  storage  of  the  WORKSET  within  the  user  program  module  is  done  by 
the  instructions  in  the  user  program  envelope  maintained  in  the  file 
ISEILERV.PHDjUSER.BAS  (see  How  to  Write  a  User  Program  Module,  p  35). 

Possible  errors  and  recovery  options; 

If  an  error  is  made  in  the  <user  file  name>,  or  no  such  file  exists  in 
the  user's  directory,  an  error  message  will  be  printed: 

ERROR  (n)  IN  LINE  ra 

error  explanation 

ABORTING  USER  PROGRAM  CALL  FOR  <u8er  file  name> 


A  return  to  the  main  MENU  will  be  made  at  this  point. 

If  a  fatal  error  occurs  during  the  execution  of  the  user  program,  then  the 
exit  to  return  to  PHD  may  not  be  made  or  not  be  made  properly.  In  such  a  case 
the  temporary  file  USER.TMP  will  continue  to  exist  on  mass  storage.  The  next 
time  that  PHD  is  initiated,  it  will  find  this  file  and  attenq>t  a  recovery, 
resulting  in  a  spurious  recovery  message  and  the  retrieval  of  a  WORKSET 
unexpectedly. 

If  for  some  reason  the  temporary  file  does  not  exist  when  the  user  program 
is  executed,  the  message  NO  DATA  BASE  TRANSFERRED  FROM  PHD — EXITING  TO  PHD  is 
given,  and  an  atteaq>t  to  return  to  PHD  is  made. 

Date  of  last  module  update;  27  Dec  82 
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BASIC  VARIABLE  MAPS  OF  PHD 

List  of  PHD  variables  accessible  froni  user  written  modules*:  VAZ/VMS 


Real 


M(9)  VARIABLE  METRICS  are  stored  in  this  array. 

X(IOOO)  Observations  are  stored  in  this  array.  The  values  for  the  ith 

observation  are  stored  in  the  same  order  as  the  VARIABLES  in  the 
BASE,  beginning  with  the  first  in  the  (i-l)*NXth  element  of  X(). 

Integer 

L%(500)  This  array  is  a  scratch  pad  for  sorting,  linking,  etc. 

N%  Contains  the  current  number  of  VARIABLES  per  observation,  i.e., 

the  number  of  VARIABLES  in  the  WORKSET. 

N9%  Contains  the  number  of  observations  in  the  WORKSET, 

String 

At(500)  Contains  the  PSN:ID  for  the  observations.  The  PSN:ID  for  the 
ith  observation  is  in  the  (i-l)th  element. 

Fi  Contains  the  WORKSET  name.  40  characters  maximum  length 

transferred  from/to  PHD. 

C$  Contains  the  current  COMMENT.  80  characters  maximum  length 

transferred  from/to  PHD. 

U$(9)  Contains  the  UNIT  names  for  the  VARIABLES  in  the  WORKSET.  16 
characters  maximum  transferred  from/to  PHD, 

U9$  Contains  the  directorv  of  the  version  of  PHD  being  used  by  the 

user  written  module  (normally  [SEILERV.PHD] ) . 

X$(9)  Contains  the  NAMEs  for  the  VARIABLES  in  the  WORKSET.  16 

characters  maximum  transferred  from/to  PHD. 

Y$  Argument  string  passed  to/from  PHD  on  exit/entry  to  the  user 

written  module  respectively.  80  character  maximum  length 
transferred  to/from  PHD. 


♦PHD  depends  on  the  values  of  these  variables  being  preserved.  They  should  be 
used  by  the  user  in  his  modules  only  with  caution. 
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List  of  vsrisbles  used  by  the  PHD  nodules*;  ISZll-M  Version 
Real: 

M(18)*  Variable  METRICS  are  stored  in  this  array. 

Z(IOOO)  Observations  are  stored  in  this  array 

Z  This  variable  is  used  as  a  pseudo  variable  to  set  the 

linewidth  of  the  terminal  in  PHD. 

Integer: 

CZ 

PZ 

IZ 

JZ 

KZ 

LZ 

LZ(500) 

MZ 

NZ*  Number  of  VARIABLES  in  current  BASE. 

N9Z*  Number  of  observations  in  WORKSET. 

QZ 

P9Z 

sz 

TZ 

String ; 

Ai 

A$(500)*‘8*  Observation  PSNtlDs  stored  in  this  virtual  array. 

*  PHD  depends  on  having  the  values  of  these  variables  preserved.  They  should 
be  used  in  a  USER  written  program  only  with  caution. 
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List  of  variables  used  by  the  PHD  modules  (Cont*d)*:  Version 

Ci*  This  string  contains  the  COMMENT  for  the  WORKSET. 

E^*  This  variable  is  used  to  store  the  ESCape  character  used  to 

control  the  terminal  functions. 

Fi*  This  string  contains  the  name  of  the  WORKSET. 

Fli 

F2$ 

F3$ 

F4i 

This  string  contains  the  option  list  used  by  the  MENU  module 
to  display  and  check  the  validity  of  option  selections. 

Tj*  This  string  contains  the  banner  (Program  for  Handling  Data) 

used  by  the  MENU. 

X$ 

Xj(18)*  VARIABLE  NAMES  are  stored  in  this  array. 

Ud(18)*  VARIABLE  UNITs  are  stored  in  this  array. 

U9$*  This  string  is  used  to  contain  the  mass  storage  device  and 

directory  where  the  PHD  module.s  are  stored.  Used  for  transfer 
to  all  modules  except  the  USER  programs. 

Yj  Argument  list  carried  from  module  to  module  in  this  -tring. 


*  I’llU  depends  on  having  the  values  of  these  variables  preserved.  They  should 
be  used  in  a  USER  written  program  only  with  caution. 
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HOW  TO  WRITE  A  USER  PROGRAM  MODULE 

Before  beginning  a  user  program  you  must  understand  the  structure  of  the 
database.  That  is,  the  BASIC  names  of  the  variables  tdiere  the  database  is 
stored  (i.e.,  the  BASE  definitions,  NAMEs,  UMITs,  METRICS,  WORKSET  NAME, 
COMMENT,  data  PSNtIDs,  and  the  data  itself).  The  table  of  PHD  variables 
accessible  from  user  written  modules  is  a  good  place  to  start.  When  the 
structure  is  understood,  the  algorithm  for  the  program  must  be  designed. 

In  the  design  of  this  algorithm  the  assumption  is  made  that  the  WORKSET  of 
PHD  is  available,  and  the  corresponding  BASIC  variable  names  must  be  used  in 
the  user  program.  In  virtually  any  user  program  it  is  necessary  to  calculate 
the  address  of  some  VARIABLE  in  an  observation.  The  algorithm  must  be 
designed  to  use  information  from  the  data  base  to  find  this  number.  For 
instance,  NZ  is  an  integer  number  that  gives  the  number  of  VARIABLES  in  the 
current  BASE. 

Using  this  value  the  user  can  calculate  the  address  of  any  particular 
VARIABLE  in  an  observation.  All  data  numeric  data  for  the  observations  are 
stored  in  the  one-dimensional  array  X().  Given  a  BASE  that  contains  ^ 
VARIABLES  per  observation,  the  element  of  X()  that  contains  the  third  VARIABLE 
in  the  12th  observation  would  be  {l2-l)*4+3-l  [using  the  formula  from  the 
table  of  variables  accessible  to  user  modules].  In  general  then  xth  VARIABLE 
of  the  yth  observation  would  be  at  (y-l)*N%+x-l .  Since  the  number  of 
observations  is  stored  in  the  BASIC  variable  N9% ,  y  could  be  checked  to  see  if 
it  lies  in  the  range  1  to  N9Z  inclusive.  The  other  elements  of  the  data  base 
that  may  be  necessary  are  the  BASE  definition  items:  NAME,  UNIT,  METRIC  of 
each  variable,  the  WORKSET  name,  COMMENT  and  the  PSN.ID  of  the  observation. 
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HOW  TO  WRITE  A  USER  PROGRAM  MODULE  (Cont'd) 

The  NAMES,  UNITs  and  METRICS  for  the  NX  VARIABLES  in  the  BASE  are  stored 
in  the  eleaents  of  the  one  disKnsional  string  arrays  X^(),u4(),  »td  the  one 
dimensional  real  array  M(),  repectively,  in  the  order  in  «fhich  they  were 
defined  (unless  SWAP  operations  have  been  done).  The  first  VARIABLE  defined 
has  its  defining  NAME,  UNIT  and  METRIC  in  element  1  of  Xt(),  UtO,  and  M(), 
tespectively.  Element  0  for  these  arrays  is  undefined.  The  WORKSET  name  is 
stored  in  the  string  variable  and  the  COMMENT  in  the  string  variable  Ci. 

The  PSN:ID  for  observation  i  is  stored  in  the  i-lth  element  of  the  string 
array  AtO.  A$()  is  a  non-dynamic  string  array  (i.e.,  all  elements  of  the 
array  are  of  a  fixed  length,  8  characters  in  this  case).  If  a  PSN:1D  string 
is  shorter  than  8  characters,  it  is  padded  on  the  right  with  blanks.  If  an 
attempt  is  made  to  store  a  string  longer  than  8  characters  into  any  element  of 
AtO,  only  the  first  (leftmost)  eight  characters  are  stored.  The  various 
string  functions  available  in  BASIC  may  be  used  to  isolate,  change  or  examine 
the  PSN  or  ID  parts  of  this  string  by  using  the  colon  which  separates  them  as 
a  marker. 

A  typical  program  might  progress  through  the  data  using  a  FOR/NEXT  loop 
and  modify  a  specified  VARIABLE.  As  an  example,  let  us  write  a  program  to 
transform  any  VARIABIdi  to  the  logfbase  e)  of  that  VARIABLE.  One  program  to  do 
this  would  be  as  follows: 

(the  unnumbered  lines  are  comments  explaining  what  the  next  line  or  lines  of 
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rao:  Version  2.2 

HOW  TO  WRITE  A  USER  PROOUVM  MODULE  (Cont'd) 

Notify  the  user  that  the  transformation  has  been  done  by  printing  a 
message  on  the  terminal. 

1080  PRINT  "TRANSFORMATION  COMPLETE" 

When  the  program  has  been  written,  a  check  must  be  made  to  make  sure  that 
the  program  lines  all  fall  in  the  range  1000-29998.  Then  the  envelope 
program,  USER. BAS,  is  added  to  the  program.  If  the  above  program  were  saved 
on  mass  storage  under  the  name  XFORM.BAS,  the  executable  user  program  file, 
XF0RM.EXE,  would  be  built  using  the  following  DCL  and  BASIC  commands  (all  user 
input  ends  with  a  □): 

^BASQ 

VAX-11  BASIC  VI. 4 
Ready 

OLD  XF0RM3 
Ready 

APPEND  (SEILERV.PHDjUSEin 
Ready 

REPLACE  XFORMD 

Ready 

EXim 

$  BAS  XFORIP 
13- INK  XFORftI 

The  executable  file,  XFORM.EXE,  would  now  exist  in  the  directory  that 
XFORM.BAS  was  stored  in.  The  program  may  be  run  by  replying  USE  XFORM  to  the 
OPTION?  request  in  the  main  MENU. 
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PARAffiTER  PILES  FOR  DATA  FITTING 

Data  fitting  in  PHD  is  done  on  the  current  WORKSBT  using  either  a  user 
witten  program  with  the  USE  option,  or  with  the  FIT  option  if  the  functional 
form  of  the  equation  to  fit  is  (or  can  be  made  to  be)  compatible  with  that 
required  by  this  option. 

The  functional  form  required  by  the  FIT  option  is  the  generalized 
polynomial  in  the  VARIABLES,  x^,  linear  in  the  parameters,  e^. 


0  -  I  (e.n  xp) 


where  the  sum  is  over  all  i  parameters  specified  by  the  parameter  file,  0. 
(i<17),  and  the  product  is  over  all  j  variables,  Xj,  in  the  WORKSET  to 
some  power  mj.  The  ln(xp  may  also  be  used  in  place  of  x^j 
A  function  is  linear  in  its  parameters  if  the  partial  derivatives  with 

respect  to  each  of  the  parameters  gives  a  set  of  equations  that  are  not 


functions  of  any  of  the  parameters.  Two  additional  constraints  are  imposed  by 
the  current  implementation  of  the  FIT  option: 

1)  The  VARIABLES  in  the  polynomial  must  be  taken  to  an  integer  power  (or 
the  natural  logarithm  may  be  used)  and, 

2)  The  exponents  must  be  in  the  range  from  -9  to  +9  inclusive. 

To  describe  the  particular  polynomial  desired,  a  file  called  a  parameter 
file  is  required  by  the  FIT  option.  The  file  may  have  any  name,  but  the 
extention  must  be  .PAR,  in  order  to  be  accessed  by  the  option.  The  file  must 
contain  the  following: 

1)  A  banner  (or  default  title  for  any  fit  to  be  made  using  this  parameter 
file)  The  banner  is  any  character  string  not  containing  a  carriage  return. 
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PARAMETER  FILES  FOR  DATA  FITTING  (Cont'd) 

2)  A  set  of  parameter  description  records,  one  for  each  parameter.  The 
parameter  description  records  contain  four  elements  separated  by  commas. 


a)  The  name  of  the  parameter.  Up  to  16  characters  (no  commas 
or  carriage  returns) 

b)  The  unit  of  the  parameter.  Up  to  16  characters  (no  commas 
or  carriage  returns) 

c)  The  initial  value  of  the  parameter.  Any  number  in  a  format 
acceptable  to  DKC  BASIC. 

d)  The  parameter  DECODE  string.  This  is  a  string  of  characters 
which  uses  positional  notation  to  define  the  function  of  each 
VARIABLE  in  the  WORKSET  to  be  used  in  the  term  containing  the 
parameter.  The  position  of  the  character  and  its  associated 
meaning  are  as  follows: 

The  first  position — Must  he  either  a  +  or  a  -.  The  + 
indicates  the  parameter  is  to  be  adjusted  by  the 
fitting  program,  ••  specifies  that  the  initial  value  of 
the  parameter  is  be  retained  during  the  fit.  If  the 
first  position  is  neither  a  nor  a  -  then  the  program 
inserts  a  -• 

All  following  positions  (2  through  the  number  of 
VARIABLES  in  the  WORKSET — The  ith  position  corresponds 
to  the  ith-1  VARIABLE  in  the  WORKSET  and  the  character 
itself  specilic's  t  lu'  function  of  the  VARIABLE  to  be 
used;  either  .1  power  of  the  VARIABLE  or  its  logarithm 
to  the  base  e. 


Ah 
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PARAMETER  FII.ES  FOR  DATA  FITTING  (Cont'd) 

The  characters  0  through  9  specify  the  zeroth  through  ninth 
power  of  the  v/ARIABLE,  respect i»^ely;  tut  letters  A-I  sj -•cify  tlife 
reciprocal  of  the  zeroth  through  ninth  power  of  the  VARIABLE, 
respectively;  an  L  specifies  the  natural  logarithm  of  the 
VARIABLE. 

For  a  valid  description  of  lIk  equation,  the-  neilier  of  positions  in  the 
DECODE  must  be  one  greater  than  the  aumber  .•!  VARlAB'JOs  in  the  WORKSET.  If  it 
is  desired  not  to  include  a  VAKIAI'I.E  in  a  tern,  the  character  0  should  be  used 
in  the  DECODE  (since  x  ”1,  this  ftectively  removes  the  VARIABLE  from  the 
term).  As  a  concrete  examr-lo,  paramote!  file  is  constructed  to  fit  the 
quadratic  equation 

2 

Resistance  =  A  +  b^Teiiperatnre  ■*  C^'Temperature 
to  a  DATASET  containing  the  variables  TEMPFRATIlRi:  and  RESISTANCE. 

Letting  R  be  resistance  and  7  be  temperature,  the  equation  is  rearranged 
to  standard  form: 

0=  A  +  B*r  +  C*'T'  +  L'i'm.  where  the  value  of  D  must  be  -1. 

If  the  DATASET  has  the  VAkiABLi's  stored  with  T  as  the  first  VARIABLE  and  R 
as  the  second  VARIABLE,  tl,.;  !)' i  .IDE,,  fvir  A,  »,  C  and  D  will  be  +00, +10, +20,  and 
-01,  respectively.  The  paramef'  r  file  must  be  created  before  fitting  using  a 
system  editor  (e.g.,  EDI  or  iiOT)  and  would  appear  as 

FIT  OF  RESISTANCE  TO  (};;AI),-AT1 IN  TEMPERATURE 

A,  OHM, 0, +00 

B, OHM/K,0,+10 

C, OHM/K2,0,+2O 

D,  OHM, -1,-01 
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•Tsr 

IV  ' 

PARAMETER  FILES  FOR  DATA  KITTING  (Cont’d) 

A 

If  the  file  in  which  these  images  are  stored  is  called  RQUAO.PAR  then  the 
fit  of  the  DATASET  name  /3T  would  be  made  using  PHD  by  GETting  the  DATASET 
RVST  and  specifying  FIT  RQUAD  in  response  to  the  next  OPTION?  request  in  the 
main  MENU. 

!  I 
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Saiq>le  PHD  USER  Program  for  the  VAX 

The  PIT  Option  written  as  a  USER  program 
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Version  2.2 


SUMMARY  SHEET  FOR  PHD  OPTION  COMMANDS 


ADD 

no  arguments  are  allowed. 

BASE 

BASE  NEH 
BASE  <data8et> 

<dataset>  is  the  name  of  a  valid  PHI)  dataset. 

CM  <optional  editing  characterXcoiiiiiient> 

<coniment>  is  any  character  string  w}iich  does  not  include  a  carriage 
return. 

Possible  <optional  editing  character>s  and  their  meanings  are 

1)  >  replace  the  rightmost  portion  of  the  COMMENT  with  this 
string.  The  replacement  is  made  on  a  character  for  character  basis. 

2)  <  replace  the  leftmost  portion  of  the  COMMENT. 

3)  -  add  this  string  preceding  the  COMMENT. 

4)  +  add  this  string  to  the  end  of  the  COMMENT. 

CORRECT 

no  arguments  are  allowed. 

DELETE 

DELETE  <sequence  seO 

<sequence  set>  is  a  set  of  <sequ(nice  range>s  separated  by  commas 

(,). 

<sequence  range>  is  <sequence  number>  or  <start  sequence 
number>-  or  <start  sequence  nuniber>-<end  sequence  number>  or 
-<end  sequence  number> 

ENTER 

no  arguments  are  allowed. 

EXIT 

EXIT  <program  file  naiK> 

<program  file  narae>  is  the  name  of  a  valid  .EXE  file  to  which  control 
is  to  be  transferred  on  exit  froi:i  P/ID. 


FIT 

FIT  <paraBeter  file  naiiie> 

<paraiiieter  file  nanie>  is  a  name  for  an  descriptor  file  that  provides 
parameter  information  for  the  fit.  The  file  extension  must  be  .PAR 

GET 

GET  <data8et  iiaBe> 

<dataset  nanie>  is  the  name  of  an  existing  dataset  on  the  mass  storage 
device  (the  file  name  is  <dataset  naraO.DAT).  The  format  of  the 
dataset  must  be  compatible  with  PUD. 
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SDIMART  SHEET  FOR  PHD  COMMAHDS  (Cont'd) 


LIST 

no  arguaents  are  allowed. 

MFP 

no  arguments  are  allowed. 

PLOT 

no  arguments  are  allowed. 

PRIHT 

PRUVT  <m^ier  of  copie8> 

■^number  of  copies>  is  an  integer  between  1  and  10. 

RENAME 

RENAME  <new  workaet  name> 

<new  workset  naine>  is  a  filename  acceptable  to  VMS. 


SAVE 

SAVE  <data8et  naae> 

<dataset  naine>  is  a  valid  VMS  file  name. 


SORT 

no  arguments  are  allowed. 

SWAP 

no  arguments  are  allowed. 

USE 

USE  ^user  file  name> 

<u8er  file  name>  must  be  the  name  of  an  .EXE  type  file  built  using  the 
USER. BAS. module  as  an  envelope.  The  user  written  program  lines  which  may 
modify  the  WORKSET  must  be  sequenced  between  1000  and  29999. 
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